Latent variable models in statistical genetics

Understanding the complexities that are involved in the genetics of multifactorial diseases is still a monumental task. In addition to environmental factors that can influence the risk of disease, there is also a number of other complicating factors. Genetic variants associated with age of disease onset may be different from those variants associated with overall risk of disease, and variants may be located in positions that are not consistent with the traditional protein coding genetic paradigm. Latent Variable Models are well suited for the analysis of genetic data. A latent variable is one that we do not directly observe, but which is believed to exist or is included for computational or analytic convenience in a model. This thesis presents a mixture of methodological developments utilising latent variables, and results from case studies in genetic epidemiology and comparative genomics. Epidemiological studies have identified a number of environmental risk factors for appendicitis, but the disease aetiology of this oft thought useless vestige remains largely a mystery. The effects of smoking on other gastrointestinal disorders are well documented, and in light of this, the thesis investigates the association between smoking and appendicitis through the use of latent variables. By utilising data from a large Australian twin study questionnaire as both cohort and case-control, evidence is found for the association between tobacco smoking and appendicitis. Twin and family studies have also found evidence for the role of heredity in the risk of appendicitis. Results from previous studies are extended here to estimate the heritability of age-at-onset and account for the eect of smoking. This thesis presents a novel approach for performing a genome-wide variance components linkage analysis on transformed residuals from a Cox regression. This method finds evidence for a dierent subset of genes responsible for variation in age at onset than those associated with overall risk of appendicitis. Motivated by increasing evidence of functional activity in regions of the genome once thought of as evolutionary graveyards, this thesis develops a generalisation to the Bayesian multiple changepoint model on aligned DNA sequences for more than two species. This sensitive technique is applied to evaluating the distributions of evolutionary rates, with the finding that they are much more complex than previously apparent. We show strong evidence for at least 9 well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least 7 classes in an alignment of four mammals, including human. A pattern of enrichment and depletion of genic regions in the profiled segments suggests they are functionally significant, and most likely consist of various functional classes. Furthermore, a method of incorporating alignment characteristics representative of function such as GC content and type of mutation into the segmentation model is developed within this thesis. Evidence of fine-structured segmental variation is presented.

European and Polynesian admixture in the Norfolk Island population

The Norfolk Island population in the South Pacific is primarily the product of recent admixture between a small number of British male and Polynesian female founders. We identified and genotyped 128 Ancestry Informative Markers (AIMs) spread across the autosomes, X/Y chromosomes and mitochondrial DNA genome, to explore and quantify the current levels of genetic admixture in the Norfolk Islanders. On the basis of autosomal AIMs, the population shows mean European and Polynesian ancestry proportions of 88 and 12%, respectively. However, there is a substantial variation between individuals ranging from total European ancestry to near total Polynesian origin. There is a strong correlation between individual genetic estimates of Polynesian ancestry and those derived from the extensive pedigree and genealogical records of Islanders. Also in line with historical accounts, there is a substantial asymmetry in the maternal and paternal origins of the Islanders with almost all Y-chromosomes of European origin whereas at least 25% of mtDNAs appear to have a Polynesian origin. Accurate knowledge of ancestry will be important in future attempts to use the Island population in admixture mapping approaches to find the genes that underlie differences in the risk to some diseases between Europeans and Polynesians.

Principal component and linkage analysis of cardiovascular risk traits in the Norfolk isolate

OBJECTIVE(S): An individual's risk of developing cardiovascular disease (CVD) is influenced by genetic factors. This study focussed on mapping genetic loci for CVD-risk traits in a unique population isolate derived from Norfolk Island. METHODS: This investigation focussed on 377 individuals descended from the population founders. Principal component analysis was used to extract orthogonal components from 11 cardiovascular risk traits. Multipoint variance component methods were used to assess genome-wide linkage using SOLAR to the derived factors. A total of 285 of the 377 related individuals were informative for linkage analysis. RESULTS: A total of 4 principal components accounting for 83% of the total variance were derived. Principal component 1 was loaded with body size indicators; principal component 2 with body size, cholesterol and triglyceride levels; principal component 3 with the blood pressures; and principal component 4 with LDL-cholesterol and total cholesterol levels. Suggestive evidence of linkage for principal component 2 (h(2) = 0.35) was observed on chromosome 5q35 (LOD = 1.85; p = 0.0008). While peak regions on chromosome 10p11.2 (LOD = 1.27; p = 0.005) and 12q13 (LOD = 1.63; p = 0.003) were observed to segregate with principal components 1 (h(2) = 0.33) and 4 (h(2) = 0.42), respectively. CONCLUSION(S): This study investigated a number of CVD risk traits in a unique isolated population. Findings support the clustering of CVD risk traits and provide interesting evidence of a region on chromosome 5q35 segregating with weight, waist circumference, HDL-c and total triglyceride levels.

Linkage disequilibrium analysis in the genetically isolated Norfolk Island population

Norfolk Island is a human genetic isolate, possessing unique population characteristics that could be utilized for complex disease gene localization. Our intention was to evaluate the extent and strength of linkage disequilibrium (LD) in the Norfolk isolate by investigating markers within Xq13.3 and the NOS2A gene encoding the inducible nitric oxide synthase. A total of six microsatellite markers spanning approximately 11 Mb were assessed on chromosome Xq13.3 in a group of 56 men from Norfolk Island. Additionally, three single nucleotide polymorphisms (SNPs) localizing to the NOS2A gene were analyzed in a subset of the complex Norfolk pedigree. With the exception of two of the marker pairs, one of which is the most distantly spaced marker, all the Xq13.3 marker pairs were found to be in significant LD indicating that LD extends up to 9.5-11.5 Mb in the Norfolk Island population. Also, all SNPs studied showed significant LD in both Norfolk Islanders and Australian Caucasians, with two of the marker pairs in complete LD in the Norfolk population only. The Norfolk Island study population possesses a unique set of characteristics including founder effect, geographical isolation, exhaustive genealogical information and phenotypic data of use to cardiovascular disease risk traits. With LD extending up to 9.5-11 Mb, the Norfolk isolate should be a powerful resource for the localization of complex disease genes.

Phenotypical characterisation of the isolated Norfolk Island population focusing on epidemiological indicators of cardiovascular disease

Objectives Only 193 people from Pitcairn Island, all descended from 9 ‘Bounty’ mutineers and 12 Tahitian women, moved to the uninhabited Norfolk Island in 1856. Our objective was to assess the population of Norfolk Island, several thousand km off the eastern coast of Australia, as a genetic isolate of potential use for cardiovascular disease (CVD) gene mapping. Methods A total of 602 participants, approximately two thirds of the island’s present adult population, were characterized for a panel of CVD risk factors. Statistical power and heritability were calculated. Results Norfolk Islander’s possess an increased prevalence of hypertension, obesity and multiple CVD risk factors when compared to outbred Caucasian populations. 64% of the study participants were descendents of the island’s original founder population. Triglycerides, cholesterol, and blood pressures all had heritabilities above 0.2. Conclusions The Norfolk land population is a potentially useful genetic isolate for gene mapping studies aimed at identifying CVD risk factor quantitative trait loci (QTL).

The B subunit of coagulation factor XIII is linked to renin and the Duffy blood group to α-spectrin on human chromosome 1

Family linkage studies were used to detect two linkage relationships on human chromosome 1. The B subunit of coagulation factor XIII showed significant linkage to renin with a maximum lod score of 5.071 at a distance of 10 cM. Significant linkage was also shown between the Duffy blood group and α-spectrin with linkage results giving a combined lod score of 3.194 at 5 cM.

Regional chromosomal assignment of human renin gene to 1q12→qter and use in linkage studies in Charcot-Marie-Tooth disease

The gene for renin, previously mapped to human chromosome 1, was further localized to 1q12 → qter using human-mouse somatic cell hybrid DNAs. The renin DNA probe used (λ HR5) could detect a HindIII restriction fragment length polymorphism. When used in studies of 12 informative families, no linkage could be found between the renin and Charcot-Marie-Tooth disease. Furthermore, an association of any renin allele with hypertension was not apparent.

Chromosome I linkage studies in Charcot-Marie-Tooth neuropathy type I

Charcot-Marie-Tooth neuropathy type 1 (CMT1) is an autosomal dominant disorder of peripheral nerve. The gene for CMT1 was originally localized to chromosome 1 by linkage to the Duffy blood group, but it has since been shown that not all CMT1 pedigrees show this linkage. We report here the results of linkage studies using five chromosome 1 markers - Duffy (Fy), antithrombin III (AT3), renin (REN), β-nerve growth factor (NGFB), and salivary amylase (AMY1) - in 16 CMT1 pedigrees. The total lod scores exclude close linkage of CMT1 to any of these markers. However, individual families show probable linkage of CMT1 to Duffy, AT3, and/or AMY1. No linkage was indicated with REN or NGFB. These results indicate that possible location of a CMT1 gene between the AMY1 and AT3 loci at p21 and q23, respectively, on chromosome 1 and support the theory that there is at least one other CMT1 gene.

Heterogeneity evidence and linkage studies on Charcot-Marie-Tooth disease

Charcot-Marie-Tooth neuropathy type 1 (CMT1) is an autosomal dominant disorder originally localized to chromosome 1 by linkage to the Duffy blood group. Studies have since shown that the disorder may be heterogeneous, as not all families show this linkage. We tested genetic heterogeneity by the HOMOG computer program in 15 CMT1 pedigrees informative for Duffy. We detected no evidence for heterogeneity in this sample, but when we combined results with previously published lod scores, heterogeneity was statistically significant. Twelve of the 15 families studied did not show linkage to Duffy. We found six of these families to be informative for a chromosome 19 marker, apolipoprotein CII(ApoC2). Despite a previous report showing probable linkage of a non-Duffy-linked CMT1 pedigree to two chromosome 19 markers, we did not detect significant linkage of ApoC2 to CMT1 in these families.

Isolation and use of chromosome 1 probes for linkage studies on Charcot-Marie-Tooth disease

Nine probes were isolated from a human chromosome 1 enriched library and mapped to regions of chromosome 1 using somatic cell hybrid lines. One clone, LR67, which mapped 1q12→q23 detected a BglI RFLP. This probe, as well as 4 other known chromosome 1 markers, α-spectrin, Factor XIIIB, DR10 and DR78, were used for linkage studies in 15 Charcot-Marie-Tooth disease (CMT1) families. Close linking of CMT1 to any of the 5 markers was not indicated. Total lod scores excluded linkage of CMT1 to LR67 and to DR10 at 5 cM or less, to DR78 and 10 cM or less, α-spectrin at 15 cM or less and Factor XIIIB at 20 cM or less. Possible linkage, however, was shown between LR67 and CMT1 at a distance of 30 cM. Also linkage at a distance of 5 cM was detected between this probe and α-spectrin.

Frequency in hypertensives of alleles for a RFLP associated with the renin gene

The genetic basis of primary hypertension is not known. Renin is important in blood pressure and volume control and a HindIII restriction fragment length polymorphism (RFLP) is present within the human renin gene locus. To examine whether there is a relationship between this RFLP and primary hypertension, DNA and renin analyses were performed on leukocytes and plasma from hypertensive and normotensive individuals. In hypertensives the frequencies of alleles for the HindIII RFLP were found to be 0.55 and 0.45, compared with 0.60 and 0.40 in the total population of 231 subjects examined, a difference that was not statistically significant. There also appeared to be no significant difference in renin activity in plasma for hypertensive patients of each genotype, nor in their pre- or post-treatment blood pressures. We thus conclude that, within the limits of the present study, the suspected genetic abnormalities associated with primary hypertension in man do not appear to be related to a HindIII RFLP in the renin gene.

DNA probes in Charcot-Marie-Tooth neuropathy

Results of Duffy (Fy) linkage confirm genetic heterogeneity in Charcot-Marie-Tooth disease type 1 (CMT1). Of 11 families informative for Fy, four showed probable linkage with CMT1, seven showed probable non-linkage and two showed definite non-linkage. These results suggest that Fy linked CMT1 may be less common than previously thought. These results combined with those of another DNA probe for the antithrombin III gene confirm that there are at least two gene loci for CMT1, termed 1A and 1B.

The cytochrome P-450 isoenzyme CYP2E1 in the biological processing of industrial chemicals : consequences for occupational and environmental medicine

The importance of the isoform CYP2E1 of the human cytochrome P-450 superfamily of enzymes for occupational and environmental medicine is derived from its unique substrate spectrum that includes a number of highly important high-production chemicals, such as aliphatic and aromatic hydrocarbons, solvents and industrial monomers (i.a. alkanes, alkenes, aromatic and halogenated hydrocarbons). Many polymorphic genes, such as CYP2E1, show considerable differences in allelic distribution between different human populations. The polymorphic nature of the human CYP2E1 gene is significant for inter-individual differences in toxicity of its substrates. Since the substrate spectrum of CYP2E1 includes many compounds of basic relevance to industrial toxicology, a rationale for metabolic interactions of different CYP2E1 substrates is provided. In-depth research into the inter-individual phenotypic differences of human CYP2E1 enzyme activities was enabled by the recognition that the 6-hydroxylation of the drug chlorzoxazone is mediated by CYP2E1. Studies on CYP2E1 phenotyping have pointed to inter-individual variations in enzyme activities. There are consistent ethnic differences in CYP2E1 enzyme expression, mostly demonstrated between European and Japanese populations, which point to a major impact of genetic factors. The most frequently studied genetic polymorphisms are the restriction fragment length polymorphisms PstI/RsaI (mutant allele: CYP2E1*5B) located in the 5′-flanking region of the gene, as well as the DraI polymorphism (mutant allele: CYP2E1*6) located in intron 6. These polymorphisms are partly related, as they form the common allele designated CYP2E1*5A. Striking inter-ethnic differences between Europeans and Asians appear with respect to the frequencies of the CYP2E1*5A allele (only approximately 5% of Europeans are heterozygous, but 37% of Asians are, whilst 6% of Asians are homozygous). Available studies indicate a wide variation in human CYP2E1 expression, which are very likely based on complex gene-environment interactions. Major inter-ethnic differences are apparent on the genotyping and the phenotyping levels. Selected cases are presented where inter-ethnic variations of CYP2E1 may provide likely explanations for unexplained findings concerning industrial chemicals that are CYP2E1 substrates. Possible consequences of differential inter-individual and inter-ethnic susceptibilities are related to individual expressions of clinical symptoms of chemical toxicity, to results of biological monitoring of exposed workers, and to the interpretation of results of epidemiological or molecular-epidemiological studies.

Carcinogenicity and genotoxicity of ethylene oxide : new aspects and recent advances

Long-term inhalation studies in rodents have presented unequivocal evidence of experimental carcinogenicity of ethylene oxide, based on the formation of malignant tumors at multiple sites. However, despite a considerable body of epidemiological data only limited evidence has been obtained of its carcinogenicity in humans. Ethylene oxide is not only an important exogenous toxicant, but it is also formed from ethylene as a biological precursor. Ethylene is a normal body constituent; its endogenous formation is evidenced by exhalation in rats and in humans. Consequently, ethylene oxide must also be regarded as a physiological compound. The most abundant DNA adduct of ethylene oxide is 7-(2-hydroxyethyl)guanine (HOEtG). Open questions are the nature and role of tissue-specific factors in ethylene oxide carcinogenesis and the physiological and quantitative role of DNA repair mechanisms. The detection of remarkable individual differences in the susceptibility of humans has promoted research into genetic factors that influence the metabolism of ethylene oxide. With this background it appears that current PBPK models for trans-species extrapolation of ethylene oxide toxicity need to be refined further. For a cancer risk assessment at low levels of DNA damage, exposure-related adducts must be discussed in relation to background DNA damage as well as to inter- and intraindividual variability. In rats, subacute ethylene oxide exposures on the order of 1 ppm (1.83 mg/m3) cause DNA adduct levels (HOEtG) of the same magnitude as produced by endogenous ethylene oxide. Based on very recent studies the endogenous background levels of HOEtG in DNA of humans are comparable to those that are produced in rodents by repetitive exogenous ethylene oxide exposures of about 10 ppm (18.3 mg/m3). Experimentally, ethylene oxide has revealed only weak mutagenic effects in vivo, which are confined to higher doses. It has been concluded that long-term human occupational exposure to low airborne concentrations to ethylene oxide, at or below current occupational exposure limits of 1 ppm (1.83 mg/m3), would not produce unacceptable increased genotoxic risks. However, critical questions remain that need further discussions relating to the coherence of animal and human data of experimental data in vitro vs. in vivo and to species-specific dynamics of DNA lesions.